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Abstract 

An attempt to incorporate relativistic kinematics in the description of light quark 
systems is made. It seems that the way of such incorporation along the suggestion 
expressed by R. Gaida and his collaborators is very promising. Comparison of these 
results with the experimental data concerning boson mass spectrum shows that this 
approach is among the best theoretical interpretation of the data. 
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Surprising success of non-relativistic quark models with funnel or oscillator potential 
in describing hadron mass spectra for heavy quark systems had waken the hopes that 
at least some aspects of the strong interactions are understood after all. Good insight 
of previous works concerning this subject is given in |l|, ^ But evidently the picture 
is not complete until it incorporates the light quarks u, d and s. But this means that 
the problem is shifting into quantum chromodynamical sector. In this respect the search 
of alternative possibilities of incorporating relativistic kinematics is quite actual. A very 
promising approach to the solution of relativistic problem of interaction of two particles 
was suggested recently by R. Gaida and his collaborators Their results can be 

directly applied to our problem of calculating the mass spectrum of quarkonium as a 
system of two-quarks. Somewhat different approach was simultaneously developed by I. 
Todorov and P. Bogolyubov ^ and later by E. Predazzi et. al. [0. In this work 
we shall investigate the results concerning the application of their findings to light-quark 
systems. 



Let us follow first the way, suggested by Predazzi et al |]T3-|112[. Following the ideas 



expressed in [|10|, |TT1 and later developed in [0, let us start with the classical expression 



for the relativistic total energy of two particle system with masses m and M, respectively 



E = ^Y>i + m? + Vpi + (1) 

were pi = — p2 = p in the cm. system. The main problem here is the mode of inclusion 
the interaction between two particles. For equal masses Lucha and Schoberl [|l^] write 
basing on (|l|) and passing to quantum mechanics 

2y^p2 + m? + V^^! (r) = Em (r) (2) 

which they call the spinless Bethe-Salpeter equation. They solve very cunningly this 
equation in the configurational momentum representation but obtain not very good re- 
sults, which we shall discuss later. But in any case we can conclude that exactly the 
way of introducing the interaction is responsible for their failure. Therefore we propose 
here another way of introducing the interaction. Namely, following Predazzi et al |T2| we 
linearize the expression 



p' + (3) 



V 

which follows from (|l]) after simple algebraic transformations, and obtain 



2E 



^ (r) = (acp + I3mc^ + V)^ (r) (4) 



where a and (3 are the usual Dirac matrices, \l/ is the four- component wave function for 
which we shall use two-component representation 



(5) 




In the general case we would have on this stage to decide what kind of Lorentz-transform 
properties we shall ascribe to the interaction. 
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In general the interaction can transform either hke Lorentz-scalar (hke the mass) or to 
be a 4-th component of Lorentz- vector, i.e. transform hke energy. Then we shah consider 
the interaction as a mixture 



V = pS ■ {l-e) + lV -e 



(6) 



where e is some mixing parameter. In what foUows we shall simplify the expression for V 
by taking e = 1/2 (which is suggested by experimental data (see,e. g. [1^ and also [p!5| ) 
what means that (IB!) can be written in the form 



V = -{(3 + l) (V + S). 



(7) 



This form is chosen in order to obtain simple non-relativist ic result V = V if there is 
no difference between S and V. Such potential was introduced previously by V. Kukulin, 
M. Moshinsky et. al. [16, l7| and was called an averaged potential, which allows to reduce 
the system of Dirac equations to a single relativistic oscillator equation. In this, particular 
case when the mixture of scalar S and vector V part of potentials are equal we shall 
consider the following possibilities 



S{r)+V{r) = V = -^ + A-r^ + Vq, 



S{r) + V{r) = V 



a 



— + k-r + Vo. 



(8) 
(9) 



S{r) + V (r) = V 



1-e 



-fir 



IGvr 



-kr 



25 



r ■ In 



+ Vo. 



(10) 



Today all of these potentials are used to describe the quark-antiquark interaction. The 
discussion concerning the advantages and handicapes of these potentials is given in 



Our aim is to apply these potentials to describing the meson mass spectra with relativistic 
kinematics which is built into equation (^. In comparison to ||18[ the search of the best 
description by minimizing will be given. For an averaged potential (^ the equation 
(H) reduces to single equation for the large wave function {p{r) 









4p2 + 4(m + ^^ 



{E - Vy ip (r) = [4p2 + 4m2 + AmS + S'^ 
For this propose we have to solve numerically the following equation 



^(r). (11) 



(12) 



Transferring to operators and carrying out the substitutions for unknown function 



(13) 



one obtains the equation 



2 



1(1 + I) [E \ V 2' 



4 



<!>(r) = 0. 



In a simple approximation 



one obtains 



S + V=UAr^ + Vo 



which leads to the equation for relativistic isotropic oscillator 



(14) 



(15) 



(16) 



dr"^ 



E \ 1 , 2 2 (E \ 1,/ 

V m] ■ -Ar H m — \- m] ■ - Vn 

2/2 4 V2 / 2 " 



Now with the standard change of variables 



E \ 1 . 2 
\- m] -A ■ r 

2 J 2 



X 



one obtains 



dx"^ 



l{l + 1) 



x^ + 



The physical solution of (0) is satisfied, when 
£;2 ' 



<l>(r) = 0. (17) 



^fr) = 0. 



_ - = . d . . + 2; + 3) + li<^/i±^ 

4 V 2 2 ^ ' 2 



or 



— m 



— - — (4iV + 2/ + 3) + — 



In the non-relativistic limit when E ~ 2m one has 



;i8) 



(19) 



(20) 



(21) 



A 

E ^2m+ \ — {AN + 21 + ?,) + Vq 
V m 



(22) 



in full accordance with the non-relativistic case (see,e.g. |T8|). 

Actually for large total energies E from (pT]) it follows E"^ ~ /^^^ i.e. almost lin- 
ear Regge-trajectory as it should be from general considerations. As we shall see the 
application of other variants of potentials gives even better results. 

As we have mentioned earlier, recently very promising vr approach to the problem 
of relativistic description of many-particle system was elaborated along with the ideas 
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presented by Gaida in |^ and elaborated by Tretyak, Spytko and Duviryak 0. They 
used the Weyl quantization method and succeeded in solving the problem for relativistic 
oscillator between two particles. 

Considering oscillator-type interaction 



V = u;^pip2r^ 



(23) 



where 



P1P2 



rriRED 



2 4 

the non-relativistic approximation they would obtain for mass of two-particle system 



(24) 



M 



2mn 



UJ 



- {AN + 21 + 3) 



2 ^ 1/2 



(25) 



or 



M 



( rx' \ A 

\ \ y2mg + J— ■ {4N + 21 + 3)] + — + Vo 



(26) 



if we express the (^) in form of string tension A and generalize their results to our 
boundary condition of isotropic oscillator. 



Exactly such approximation was taken in |T8[. But if we suggest that in general 



P1P2 



E + m„ 



then Lvov group results 



M 



2mq + 



2A 



E 



{AN + 21 + 3) 



+ 



2A 



1/2 



E + ma 



(27) 



(28) 



will resemble very much our relativistic approach. The results of calculations according 
to ( PBD are given below. The parameters are taken to be A = 0.01 GeV^, Vq = —0.436 
GeV certainly this approximation is valid if one considers the energy dependence of pj 
only on the last stage of calculations. 

It is interesting that similar to (^) result for M^- operator was obtained in relativistic 
approach by Ishida-Oda based on special assumption of covariant relativistic approach. 
Even the numerical values of parameters of Ishida-Oda are close to the result of (p8|). 
Their value A = 0.05GeV^ which is of the same order as our A. More precise comparison 
is impossible because of ambiguity of other parameters. 

Usually such potential models like we have used here are called naive quark model. But 
exactly our model is not-so naive. Firstly, the relativistic kinematics not only renders it 
more complicated, but shows the possible way of building the model of interaction of two 
relativistic particles. Secondly the potentials like (p|-p!OD incorporate asymptotic freedom 
the strong coupling constant ag was calculated according to the classical expression 



a. r 



33 - 2N 



Vr2A2 



(29) 
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where A was taken to be equal to A = O.lAGeV. And lastly the model allows to include 
spin-spin interaction either by passing to Breit-Fermi equation or by using Dirac equation 
straightly. In the Table 1 is taken exactly according to (]29|). 

It is interesting to note that the definition of mass can be given in a different way. 
Considering one particle as moving in the field of another and vice versa and adding the 
obtained masses we obtain the results with parameters which correspond very closely to 
our Table 1. Namely one can write down the Dirac equation for one particle moving in 
the outer field, reduce it to the equation for large component x (r) and obtain ^0 



(e^ - m") X (r) = p'x (r) + [E + m) Vx (r) . 
Applying the virial theorem to this equation 



x{y) v'^{y) x{y) 



x{y) 



dv 

' dy 



x{y) 



x{y) x{y) 



Combining (30) and (31) for oscillator interaction one obtains 



m 



where 



B 



Nl 



2JC {E + m)Bm + {E + m) Vo 



(x(y)|2/'|x(y)) = 2Ar + Z + 3/2 



(30) 



(31) 



(32) 



(33) 



This expression on resembles very much the (|2^). So we are left with the five most 
realistic from our point view possibilities, namely, of calculating masses according to (|]), 

m 



i), (0), (H), dH). 

To obtain the masses of multi-quark system according to these expressions, it is neces- 
sary to define the values of parameters. The masses of quarks were taken to be as usually 
in quark models = 0.33(je\^, rric = 1.675 — l.TSGeT^, rrib = 5.05 — S.lGeV. Parameters 
for different potentials are shown in Table 1. 



Table 1. The parameter values for different potentials 



Potential 






OLh 


k, A 


Vo, GeV 


i) 


0.5 


0.325 


0.3 


k = 0.27 

(GeV^) 


-0.8356 


m 


0.5 


0.386 


0.3 


A = 0.01 

(GeV^) 


(-0.436)... 
...(-0.527) 


m 


£- = 0.3795 

(GeV^) 


/i = 0.054 

(GeV) 


K = 0.75 

(GeV) 


A = 0.35 

(GeV) 
h = A 


-1.103 
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Table 2. Variation of for different parameters of potential (j^ 



k 

(GeV^) 


A,uu 




^bb 


Vo 
(GeV) 


0.1826 


514 


1.15 ■ 10^ 


1.3 ■ 10* 


-0.628 


0.25 


98 


6.2 ■ 10^ 


2 • 10* 


-0.791 


0.26 


67 


2.2 ■ 10* 


1.6 ■ 10* 


-0.813 


0.27 


43.6 


3-10* 


2.9 ■ 10^ 


-0.835 


0.29 


12.6 


7-10* 


8.5 ■ 10^ 


-0.878 


0.305 


1.8 


1.36 ■ 10^ 


4.3 ■ 10* 


-0.908 



The value of Vq reflects the fact that potential is the Fourier-transform of scattering 
amplitude Vq being the constant of interaction. Gromes [21] has evaluated this constant 
for linear confinement and obtained the value of Vq 



Vn 



-2v^-e-(^~°-5) 



(34) 



where k is the string tension 7 = 0.57721... is Euler-MacLoraint constant. According to 
our values of k = 0.18 — 0.3056*6^^. According to (0) has to vary within the limits 
Vo = — (0.77 ^ 1.02) Gel^ which is quite close to values, cited in the table (H). Let 
us stress that the set values are not colour-dependent, which, reduces the number of 
adjustable parameters. 

The result of calculation together with experimental data are shown in Tables (0-§|). 
Experimental values were taken from p2|. For choosing the parameters the minimum of 



X — criterion was used, with the definition of x given in |T8[. In this definition N is the 



number of meson masses, n— is the number of parameters (in our case we considered them 
to be equal to two namely confinement parameter and Vq), A is the experimental error 
in definition of experimental mass Mexp of two-quark system since Vq was chosen to 
match the experimental value of ground-state mass we are actually left only which one 
adjustable parameter A (or k). Since we do not include LS-forces we had to take the 
average center of gravity (COG) value of P-resonances, which was calculated according 
to formula 



COG 



Ej (2i + 1) ■ Mj 
E, (2j + 1) 



(35) 



As one can see from both radial and orbital excitation calculation the variant which 
incorporates in one or another way the relativistic kinematics give better description of 
Regge-trajectories, which, are believed to be linear in / for M^. 

It is well known that even in non-relativistic limit one can obtain good description of 
uu-dd -systems at the account of spoiling the cc of bb -description. 

Indeed in the Table 2 we show obtained for Eichten parameters. Due to a high 
precision of defining J/^> and T-mesons the large x^ were considered good for these 



mesons and bad for p-meson trajectory. On other hand the non-relativistic Badalyan p3| 
results are good for uu, but bad for J/^. Fabre |^ in order to obtain good results for 
light quarkonium had to change the very potential. Instead of this situation our way or 
incorporating the relativistic kinematics we obtain good results for all data. The table 
2 demonstrates this statement. Table 3 contains the comparison of different potentials. 
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The discussion concerning the choice of potentials is given in We shall choose in 

what follows the Cornell-potential (P) which seems to be preferable though potential ([T0|) 
is also quite good. The parameters here are taken from table 1. The uu-data are 
fantastically good, but cc and 66-data could be better. Therefore in tables 4-6 we give the 
results for k = 0.296*6^^. We consider these results as the best. It is interesting that the 
values which give these best results are close to those of Lucha and Schoberl .We want 
to stress that all the above results are obtained by numerical solution of ([I^). The last 
columns in (4 — 6) are calculated according to (|2^ ). With the choice A = 0.071GeV^, 
Vo = — 1.0077Ge\^, niu = 0.33GeV^, rric = l.75GeV, rrif, = 5.136*61^ the results are quite 
comparable with the other entries. But still we have to conclude that pure oscillator 
potential is too rough to give final result. The use of more sophisticated potential is to 
be taken hear too. But it demonstrates nicely that inclusion of relativistic kinematics is 
very crucial. 

We would like to call the attention to one interesting feature of relativistic models, 
namely that the slope of linear (or close to linear) Regge-trajectory in this case is constant, 
while in nonlinear models it is neither constant nor linear. Roughly experiments show 
this slope to be equal ~ 1.2GeV^. In our cases it varies from 1.15 to 2.5GeV^ for different 
cases. As Tutik et ||2^ have indicated the Regge trajectories for low-lying states coincide 
while for large values of orbital momentum / the screened potential ([T0|) leads a limited 
Regge-trajectory in contrast to infinitely rising trajectories for other potentials like (W) or 



Table 3. Mass spectrum of uw-system with some realistic potentials. 



State 


designation 


Mexp 
GeV 


Mth® 
GeV 


Mth® 
GeV 


Mth& 
GeV 


IS 




0.768 ± 0.0005 


0.768 


0.768 


0.768 


IP 


3 p* 
^COG 


1.262 ±0.03 


1.26 


1.296 


1.3 


2S 




1.465 ±0.025 


1.609 


1.573 


1.561 


2P 


a2 l-(2++) 


1.935 ±0.015 


2.003 


1.932 


1.889 


3S 


pl+(l-) 


2.15 ±0.01 


2.3 


2.155 


2.103 


ID 


P3 l+(3-) 


1.691 ±0.013 


1.665 


1.689 


1.687 


2D 


P3 l+(3~) 


2.25 ±0.01 


2.348 


2.235 


2.1845 


IF 


04 l-(4++) 


2.037 ±0.036 


2.03 


2.021 


2.003 


IG 


P5 l+(5-) 


2.350 ±0.015 


2.366 


2.312 


2.275 


IH 


as 1"(6++) 


2.45 ±0.13 


2.685 


2.576 


2.514 


x' 






53.3 


1.8 


19.6 
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Table 4. The best mass spectrum of -u-u-system with Cornell-potential (Qj 
in formula p^). 



Sto,tp 


dp sionofAon, 


Mexp 
GeV 


GeV 


Mth^ 
GeV 




pl+(l--) 


0.768 ± 0.0005 


0.768 


0.768 


IP 


^COG 


1.262 ±0.03 


1.281 


1.198 


25 


p 1 + (1") 


1.465 ±0.025 


1.551 


1.577 


2P 


a2l-(2++) 


1.935 ±0.015 


1.902 


1.923 


35 


pl+(l-) 


2.15 ±0.01 


2.12 


2.244 


ID 


P3 l+(3-) 


1.691 ±0.013 


1.665 


1.577 


2D 


PS l+(3--) 


2.25 ±0.01 


2.197 


2.244 


IF 


04 l-(4++) 


2.037 ±0.036 


1.987 


1.923 


IG 


P5 l+(5--) 


2.350 ±0.015 


2.278 


2.244 


IH 


as l-(6++) 


2.45 ±0.13 


2.529 


2.547 








12.6 


35.68 



Table 5. The same for cc-system. 



State 


designation 


Mexp 
GeV 


Mth® 
GeV 


MTHm 
GeV 


IS 


j/^ o-(i~) 


3.096 ± 0.00009 


3.07 


3.014 


IP 


Xci o+(i++) 


3.51 ±0.00012 


3.5118 


3.331 


ID 


^ ?-(l— ) 


3.770 ±0.0025 


3.837 


3.632 


25 




3.688 ±0.0001 


3.732 


3.632 


2D 


* ?-(l— ) 


4.159 ±0.02 


4.296 


4.195 


35 


* ?-(l~-) 


4.04 ±0.01 


4.227 


4.195 


3P> 


m v\i—) 


4.415 ±0.006 


4.692 


4.717 








7-10^ 


8.4-105 



Table 6. The same for 66-system. 



State 


designation 


Mexp 
GeV 


Mth® 
GeV 


Mth(28) 
GeV 


15 


T ? ■• (1—) 


9.460 ± 0.00022 


9.479 


9.542 


IP 


3 p* 
^COG 


9.892 ±0.0007 


9.883 


9.741 


25 


T ?'(1"") 


10.023 ±0.00031 


10.037 


9.932 


2P 


3 p* 
^COG 


10.268 ±0.00057 


10.299 


10.121 


35 


T ?'-(l--) 


10.355 ±0.0005 


10.433 


10.308 


45 


T ?-(l"") 


10.58 ±0.0035 


10.776 


10.675 


55 


T ?'-(l--) 


10.865 ±0.008 


11.062 


11.034 


65 


T ?'(1--) 


11.019 ±0.008 


11.332 


11.385 


x' 






8.5 ■ 10^ 


7.4 ■ 10^ 
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